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Abstract. The NJL-jet model provides a sound framework for calculating the fragmentation func- 
tions in an effective chiral quark theory, where the momentum and isospin sum rules are satisfied 
without the introduction of ad hoc parameters [1]. Earlier studies of the pion fragmentation func- 
tions using the Nambu-Jona-Lasinio (NJL) model within this framework showed good qualitative 
agreement with the empirical parameterizations. Here we extend the NJL-jet model by including the 
strange quark. The corrections to the pion fragmentation function and corresponding kaon fragmen- 
tation functions are calculated using the elementary quark to quark-meson fragmentation functions 
from NJL. The results for the kaon fragmentation function exhibit a qualitative agreement with the 
empirical parameterizations, while the unfavored strange quark fragmentation to pions is shown to 
be of the same order of magnitude as the unfavored light quark's. The results of these studies are 
expected to provide important guidance for the analysis of a large variety of semi-inclusive data. 
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1. QUARK FRAGMENTATION AND NJL-JET MODEL 

Quark fragmentation functions have long been of interest for analyzing hard scattering 
reactions [2, 3, 4, 5, 6, 7, 8]. New experimental efforts for extraction of fragmentation 
functions from deep-inelastic lepton-nucleon, proton-proton scattering and e + e~ anni- 
hilation data [9, 10] have generated a renewed interest in the subject. The analysis of 
the transversity quark distribution functions [7, 11] and a variety of other semi-inclusive 
processes [12, 13] also critically depend on the knowledge of the fragmentation func- 
tions. 

The NJL-jet model of Ref. [1] provides a self-consistent framework for calculations of 
both quark distributions and fragmentation functions in an effective chiral quark theory. 
The advantage of the model is that there are no ad hoc parameters introduced to describe 
fragmentation functions; the quark self-energy normalization factor in the coupled inte- 
gral equations for the fragmentation functions arises naturally from the product ansatz. 
In this work we extend the NJL-jet model by introducing the strange quark, thus allow- 
ing fragmentation to both pions and kaons. The inclusion of the new channel is shown 
to bring significant corrections to the previously calculated fragmentation functions to 
pions. 



2. INCLUDING THE STRANGE QUARK 



2.1. Strange Quark Constituent Mass and Coupling 

Introducing the strange quark in the NJL-jet model involves calculating the strange 
quark distribution and fragmentation functions in the NJL framework, requiring the 
knowledge of the quark-meson (strange-light to kaon) coupling constant and the strange 
quark constituent mass. We calculate the quark-meson coupling using the same approach 
as Ref. [14]. The strange quark mass is chosen semi-empirically to best fit the model 
calculated kaon decay constant to the experimental value. 

The quark-meson coupling constant is determined from the pole in the quark- 
antiquark t-matrix at the considered meson's mass. This involves the derivative of the 
familiar quark-bubble graph: 
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Evaluating the integral over q + using the complex residue theorem and making the 
variable substitution q- = xk- yields for the quark-meson coupling: 
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In this work we use the Lepage-Brodsky (LB) "invariant mass" cut-off regularization 
as described in Refs [14, 1]. Using the experimental value of fx = 0.114 GeV yields 
a strange constituent quark mass M s = 0.45 GeV and the corresponding quark-kaon 
coupling constant of gK qq = 3. 



2.2. Quark Distribution and Fragmentation Functions 

The quark distribution function f q (x) has an interpretation as the probability to find 
a quark of type q with momentum fraction x in the hadron (in our case meson) h. The 
corresponding cut diagram is shown in Fig. la, which can be equivalently represented 
by the Feynman diagram depicted in Fig. lb: 
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FIGURE 1. Quark distribution functions. 
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FIGURE 2. Quark fragmentation functions. 



where k_ = xp- and Cj is the corresponding flavor factor (Cj = 2 for all the mesons 
considered except for Kq, where C/ = 1). 

The elementary fragmentation function depicted in Fig. 2 can be written as: 



W = N c ^gl qq y^Tr[S l (k)Y + S 1 (k)r 5 (li-^+M 2 )r 5 ] 
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3. GENERALIZED NJL- JET 

The NJL- Jet model of Ref. [1] uses a multiplicative ansatz for the total fragmentation 
function to derive an integral equation for the quark cascade for the process depicted in 
Fig. 3. The derived integral equation for the total fragmentation function is: 
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FIGURE 3. Quark cascade. 
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Here <i™(z) = d™(z)/(l — Zq), where Zq is the residue of the quark propagator in 

the presence of meson cloud (see Ref. [1]). Then £ m / d™(z)dz = 1, thus allowing an 
interpretation as the probability of an elementary process. In Eq. (7) the sum is over 
the flavor of the emitted quark Q and the splitting function of quark q into Q with a 
momentum fraction z is naturally the same as the splitting function of q into a meson m 
of a flavor composition qQ with a momentum fraction 1 — z. 

The result in Eq.(7) resembles the integral equation ansatz of Field and Feynman's 
quark-jet model [15, 2]. Rewriting the above expression helps to elucidate the proba- 
bilistic interpretation of the model: 

D™(z)dz = <%(z)dz + Y, f l dQ(y)dyir Q ( Z -) d ^. (8) 
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Here the left hand side term has the meaning of the probability to create a meson 
m carrying the momentum fraction z to z + dz of initial quark q. The first term on the 
right hand side corresponds to the probability of creating the meson with momentum 
fraction z to z + dz in the first step of the cascade, plus the second term corresponding 
to the creation of the meson further down the quark cascade after a splitting to a quark 
Q with momentum fraction y. Here the probability of creating the meson m scales with 
the momentum fraction left to the quark after the splitting z/y, which is clearly only the 
case in the Bjorken limit. Thus the model can be trivially generalized by including the 
strange quark directly in the Eq. (7). 

We solve the coupled set of integral equations using the elementary fragmentation 
functions of Eq. (7) for u, d and s quark fragmentation to a given meson. The cor- 
responding fragmentation functions for the anti-quarks are obtained using the charge 
symmetry from fragmentation functions of quarks to the corresponding anti-meson. The 
comparisons with the phenomenological parametrizations of Ref [9] are performed by 
DGLAP evolving the calculated fragmentation functions from the low-energy model 
scale of Q 2 = 0.18 GeV 2 to 4 GeV 2 at leading order using the software from Ref. [16]. 
The evolution kernel for the distribution functions is modified as described in Appendix 
B of Ref. [1] 

It is easy to see, using the properties of d™ along with the normalization condition of 

d™, that the solutions of the integral equations (7) should satisfy both momentum and 
isospin sum rules. Our numerical solutions obey these rules within numerical errors of 
less than a percent. 

4. RESULTS AND CONCLUSIONS 

The results for the fragmentation functions of u, d and s quarks to iz + and K + at the 
model scale are shown in Fig. 4. Though the fragmentation functions of the unfavored 




strange and light quarks are of the same order of magnitude, we can see a notable 
difference between them, even at the low scale of the model. Thus one introduces 
considerable and ultimately unnecessary uncertainties by simply assuming that ZH = 

Df , etc. (for example as is done in Ref. [16]) . 

Further we present our results for the fragmentation functions Df { and Df t = Df { in 
Fig. 5. Here the DGLAP evolved curve is also compared to the empirical parametriza- 
tions of the experimental data of Ref. [9], evolved to the same scale. We can see that 
the inclusion of strangeness softens the high z region of Df t compared to previous cal- 
culations [1], thus bringing the curves closer to the phenomenological parametrizations. 
This is expected as the elementary fragmentation to a Kaon is not negligible. In fact the 
fragmentation to K + is about half as likely as that to n + as can be seen from the plots 
in Fig. 6. The plots show a reasonably good agreement with the parametrizations within 
the large uncertainties of the latter. 

It is clear that for a more complete description of the quark fragmentation both vector 
meson and nucleon anti-nucleon channels need to be included in the calculations. This 
can be accomplished within the current framework. The high z region of fragmentation 
functions are dominated by "few-step" transitions where the availability of the additional 
fragmentation channels might have a noticeable effect. 

Another limitation of the model is the assumption of the momentum scaling of the 
probability of hadron creation in each step of the decay chain, which is clearly only 
the case in the Bjorken limit. A quark with a finite momentum loses energy with each 
production of a hadron and finally recombines with the remnants of the antiquark jet to 
form the final hadron. A more accurate description of the process requires Monte-Carlo 
(MC) simulations of the quark fragmentation, similar to the studies in Refs. [15, 17], 
and others. MC simulations would also allow one to access the transverse momentum 
distribution of the produced hadrons, thus becoming relevant for the analysis of a large 
variety of semi-inclusive data. 
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FIGURE 6. Kaon fragmentation functions. 
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